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Single electron pumping based on the topological property of Majorana fermions (MFs) is pro- 
posed. The setup consists of a quantum dot and four nano topological superconductors (TSs) 
connected by constriction junctions, with an additional vortex located in the loop of TSs. Op- 
eration is performed by gate voltages at constriction junctions. Simulations with Bogloliubov-de 
Gennes equation demonstrate successfully quantum protection during switching operation. 
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I. INTRODUCTION 

In the electronics-based computers, computation is 
performed based on motions of thousands of electrons, 
where the accuracy of bit manipulation relies on the 
statistics of huge number of electrons. If one can control 
electrons one by one at given instants, a single electron 
bit may be achieved, which reduces the energy consump- 
tion to the limit. This technology is also important for 
forming the metrology triangle of Ohm's law: the quan- 
tum Hall effect discovered by von Klitzing works as the 
conductance standard[l|, and the Josephson effect pro- 
vides the voltage standard^, while realization of precise 
current standard remains to be a challenge [3j- 

Manipulation of individual electrons in condensed mat- 
ters is not easy since the wave length of electrons is com- 
parable with their separation. One way to transport sin- 
gle electrons is to use quantum tunneling phenomenon. 
One needs to puncture the time for electrons to tunnel 
through the barrier by introducing several quantum dots 
(QDs) in the circuit, since quantum tunnelings normally 
take place randomly in time [3, [!]■ Several schemes have 
been proposed for single electron pumping so far, such as 
single electron transfer in metallic nanostructures based 
on the combination of Coulomb interaction between elec- 
trons and their quantum tunneling trough an insulating 
barrier Q, adiabatical charge pumping through a quan- 
tum dot by driving two independent parameters of the 
system coherently, e.g. phases of gate voltage and mag- 
netic field @, IH, quantum Hall effect in Corbino disk ge- 
ometry combined with a time dependent flux |9j. 

In a topological superconductor characterized by zero- 
energy Majorana fermions (MFs) as quasiparticle excita- 
tions, the ground state exhibits degeneracy with respect 
to the parity of electron number [TOl - 1261 ]. This property 
provides a new principle for single electron pumping with 
unprecedented precision guaranteed by topology. In the 
present work we formulate a protocol for this purpose. 

Our setup consists of four topological superconductors 
(TSs) each carrying a vortex and connected by constric- 
tion junctions; there is an additional vortex inside the 
loop of TSs (see Fig. 1). In order to reveal the topo- 




FIG. 1: Schematic device setup of a NOT gate for MF qubit. 
Finite square samples (called bricks) of heterostructure of fer- 
romagnetic insulator and spin-orbit coupling semiconductor 
are positioned on the surface of an s-wave superconductor 
with vortices (black dots) trapped in it. The four bricks (on 
the right side) form the NOT gate. The edge MF of another 
brick (on the left side) can be driven to circle the vortex at 
the center of the loop formed by the four bricks by tuning 
gate voltages at the point-like constriction junctions. 



logical property of the system, we first attach another 
TS to the device. We show that the edge MF of this 
TS can be driven in a controlled way around the loop 
of the four TSs by switching on and off gate voltages at 
constriction junctions in the designed sequence. After 
this process, the parity of this TS is nipped since the 
edge MF acquires a 7r phase due to the vortex at the 
center of device. Therefore, one can consider this TS 
as a MF qubit, and the device of four TSs as a NOT 
gate. We then demonstrate that the NOT gate works 
for single electron pumping when a QD with Coulomb 
blockade effect is attached when the energy of QD is 
adjusted appropriately. Numerical simulations based 
on the time-dependent Bogloliubov-de-Gennes (TDBdG) 
equation 27] are performed, which confirms successfully 
the quantum protection and phase coherence during the 
whole process typically of several nano seconds. 
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The remaining part of this paper is organized as fol- 
lows. We discuss in Sec. II the dynamics of an edge MF in 
a MF-qubit driven through the NOT gate. Then we re- 
veal the relation satisfied by the interactions in the NOT 
gate in Sec. III. Based on the property of NOT gate, 
we formulate in Sec. IV the topological single electron 
pumping by attaching a QD to the NOT gate, and give 
explicitly the energy regime for QD. In Sec. V, it is clari- 
fied that the function of NOT gate can be understood as a 
quantum interference between two MFs. Discussions are 
presented in Sec. VI with a summary given in Sec. VII. 

II. NOT GATE FOR MF-QUBIT 

A topological superconducting state can be achieved 
in a heterostructure of s-wave superconductor (s-SC), 
spin-orbit coupling semiconductor (SOSM), and ferro- 
magnetic insulator (FMI) with one vortex at the sample 
center, where one MF appears at the vortex core and 
another MF at the sample edge [22|, H3] • In the present 
device, four finite TSs are positioned on a common s-SC 
substrate. The core MFs are stable and do not partici- 
pate in the phenomena discussed below (strictly speaking 
there are exponentially small contributions, which can be 
neglected safely), and thus will be omitted hereafter. The 
linear dimension of TSs should be in the regime of tens 
to one hundred nano meters[27]: for too large TSs, the 
energy gap between the zero-energy MF state and the 
lowest excited state becomes very small, which limits the 
operation temperature, while for too small TSs, core MF 
and edge MF interact with each other, which destroys 
the MF ground state. Edge MFs interact with neighbor- 
ing ones through the constriction junctions when electron 
hoppings are permitted. There is an additional vortex in 
the s-SC substrate inside the loop formed by the four 
TSs, which governs the interactions among edge MFs as 
will be revealed below. We attach a MF qubit to the 
above device as schematically shown in Fig. 1. 

The low-energy physics of the system can be described 
by the following effective Hamiltonian, 

Hmf = iA (t)r 7o7i + i Y ^'(*) r nj'7i+i (1) 

3 = 1,4: 

where 70 denotes the edge MF at the qubit, and 7, 
with 1 < j < 4 denote those on the TSs in the NOT 
gate and 75 =71; time-dependent, dimensionless factors 
< \j(t) < 1 are introduced for the switching process 
of constriction junctions, whereas Tj are the MF interac- 
tions when the gate voltage is off and thus two TSs are 
connected fully. 

At the initial stage, the switch configuration is de- 
scribed by the vector (Ao, Ai, A 2 , A3, A 4 ) = (0,1,0,1,0) 
(see the left inset of Fig. 2). Since the qubit is isolated, 
there is an edge MF (70) localized at the qubit. In con- 
trary, with Ai = A3 = 1, both the unified edge of TS(1) 
and TS(2) and that of TS(3) and TS(4) contain two vor- 
tices, and thus there is no edge MF in the NOT gate. In 



the representation of [1] the MFs 7, for 1 < j < 4 are 
fused to finite energies due to the interactions. 

Turning on the connection between the qubit and 
TS(1), namely Ao = — > 1, the wavefunction of 70 
spreads to the unified edge of the now connected qubit, 
TS(1) and TS(2), since it contains three vortices. We 
then turn off the connection between TS1 and TS2, 
namely Ai = 1 — ► 0. The wavefunction of edge MF col- 
lapses totally on TS2 due to the topological property. 
After these two switchings, the edge MF 70 is trans- 
ported completely to TS2. Repeating this process, one 
can drive the edge MF 70 through the NOT gate in a 
clockwise way, and returns it back to the initial position 
at the qubit, with the switching sequence (0, 1, 0, 1, 0) h-> 
(1,1,0,1,0) ^ (1,0,0,1,0) h+ (1,0,1,1,0) 1 ^ 
(1,0,1,0,0) 1 — ^ (1,0,1,0,1) 1 y (0,0,1,0,1). During this 
process the edge MF feels the gauge field formed by the 
central vortex, and thus acquires a phase of n which 
makes 70 i-> —70- As the result, the electronic parity 
of the qubit is flipped. 

In order to confirm the function of the NOT gate, we 
perform numerical calculations based on TDBdG equa- 
tion. We first diagonalize the tight-binding BdG hamil- 
tonian [27j of the system in the initial stage, 

Ho = -t ^2 a L a J^ - v Y + v *(4t £ n- - a l a u) 

ij,cr i,cr i 

+ it a ~ c\ + 8 y °xCi + h-c^ + ^ [a^c^ + h.c 

i.S i 

(2) 

where both spin-conserved hopping to and spin-flipped 
hopping t a are between nearest neighbors with a the 
grid spacing, and V z is the Zeemann energy. An edge 
MF state is obtained at the MF qubit with the wave 
function \^(t = 0)) = |</>mf)- Then we modulate dynam- 
ically the hopping parameters at the constriction junc- 
tions, which changes the MF interactions in [1] sequen- 
tially and drives the edge MF % [13] ■ The evolution 
of the wave function is obtained by solving the TDBdG 
equation ih^\^(t)) = H\^(t)) based on the Chebyshev 
polynomials expansion 29, 30(. To monitor the evolution 
of the edge MF, we project the wave function onto the ini- 
tial one, and evaluate the parameter 0(t) = (</>mf| , F(£))- 
As can be seen from Fig. 2, 0(t) changes from positive 
unity at the initial stage to negative unity at the final 
stage, representing a sign change in the MF wave func- 
tion. It is worth noticing that the conservation of the 
function norm as seen in Fig. 2 confirms the topological 
protection of the edge MF during the operating process. 

III. PARITY FLIPPING IN NOT GATE 

Since the parity of the whole system is conserved upon 
application of gate voltage as well as Cooper pair tun- 
neling from SC substrate, the above switching operation 
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FIG. 2: (Color on line) Time evolution of the MF wave func- 
tion transported through the NOT gate in terms of projec- 
tions 0(t) to its initial wave function. The time for one 
switching step is T = 4 x 10 4 ft/to- The two insets schemati- 
cally show the initial and final stages, where the solid (dotted) 
lines denote the on, X — 1 (off, A = 0) state of the constriction 
junctions between TSs. 



should reverse the parities of the MF qubit and NOT 
gate simultaneously. In order to check the electronic par- 
ity of the NOT gate, we investigate first the signs of MF 
interactions 1% denned in Eq. (1). As revealed in the 
previous work[23|, when 70 is transported to TS(2), it 
picks a sign sgn(FoFi). Therefore, the sign of edge MF 
70 after driven through the NOT gate is given in terms 
of the interactions by 



7o 



sgn[r rir2r 3 r4(-ro)]7o. 



(3) 



The minus sign attached to Tq is due to the opposite 
motion of MF against the direction used to define the 
interaction To in [lj27j . 

It is then clear that the sign reversal of 70 implies 
Iil^rsI^ > 0, a topological property generated by the 
central vortex. Hereafter, we consider explicitly the case 
where all interactions are positive since all possible con- 
figurations of interaction signs can be transformed to 
each other by gauge transformation. It is worth noting 
that the same sign constraint rij=i 2N > and gauge 
choice are available for even number of TSs, which will 
be used for discussions below. 

We define two regular electronic states with the four 
MFs, d\ — (71 + «72)/2 and d\ = (73 + i^)/2, as always 
possible even when the MFs are bounded and not free. 
We then rewrite the Hamiltonian [T] in terms of the ba- 
sis \mn 2 ) = {|00), |11), 1 10) , |01)}, with m denoting the 
parity of the electronic state, 



Hmf = 



r 3 r 2 -r 4 

-(n+n; 



— v 
2 1 4 



ri-n n, + r 4 
n, + r 4 -(r;-r 3 ) 



where Tj = Xj(t)Tj. The four eigen energies are given 

by e 1>2 = ±V(ri + r 3 )2 + (r^ 
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FIG. 3: (a) Topological single electron pumping realized by 
the NOT gate and QD. A gate voltage (yellow cylinder) is 
used to tune the coupling between the QD and one TS of 
the NOT gate, (b) Working mechanism of single electron 
pumping (see text). The two lowest energy levels (red and 
blue) associated with opposite parities of the NOT gate cross 
each other upon the switching operation, with AE the energy 
difference at the final stage. When AE > en, the occupation 
energy of QD, an electron is emitted from the NOT gate to 
the QD to lower the energy of whole system. 



r 4 ) 2 for even parity, 



and E 3A = ± y / (T[ - T' 3 ) 2 + (T' 2 + T 4 ) 2 for odd parity 
The switch configuration of the constriction junctions in 
the NOT gate changes after the operation: Ai = A3 = 1 
and A2 = A4 = at the initial stage, while Ai = A3 = 
and A2 = A4 = 1 at the final stage. It is clear that 
the ground-state energy E s = —(Ti + F3) at the ini- 
tial stage is achieved in the even-parity subspace, while 
Eg = — (T2 + T 4 ) at the final stage in the odd-parity 
subspace. Therefore, the electronic parity of the NOT 
gate is reversed after the operation of switching. It is 
easy to see that this parity reversal is realized because 

r 1 r 2 r 3 r 4 > o. 



IV. SINGLE ELECTRON PUMPING 

The NOT gate discussed above can be used for single 
electron pumping when a QD is attached (see Fig. 3(a)). 
The parity flipping specifies an odd number of elec- 
trons transferring between the NOT gate and QD upon 
the switching operation described above, whereas the 
Coulomb blockade effect on the QD limits the charge 
transfer to a single electron. 

Let us first consider the case of a vacant QD with an 
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FIG. 4: (Color on line) Interference of two MFs in a loop of 
2N TSs with a central vortex. First two MFs are created at 
TS(1) and TS(2) by turning off Ai, which is denoted by the 
cross. These two MFs are then transported to TS(21) and 
TS(21+1) by turning on and off the constriction junctions on 
the paths. Finally, they are fused by turning on X21 . 

energy level eo in the initial stage (see Fig. 4(b)). The 
ground state of the total system exhibits even parity, 
since the NOT gate is in a state of even parity. From 
the parity conservation of the total system, there are two 
candidates for the ground state at the final stage after 
switching operation discussed above: vacant QD with 
total energy — |r 2 — r 4 |, and/or occupied QD with total 
energy — (T2 +T4) + €q, both exhibit even parity. One 
electron on the NOT gate is transferred to the QD upon 
the switching operation if — (r 2 + T4) + eq < — |ra — T^J, 
namely eo < 2min{r2,r4}. Physically the inequality 
means that if the occupation energy of the QD is too 
large, electron will not jump to QD. In the same way, 
one can figure out that when QD is occupied at the initial 
stage, an electron is transferred to the NOT gate upon 
switching operation if — (r 2 + T 4 ) < — |T 2 — T 4 | + e , 
namely eo > — 2 min{r2, The physical meaning of 
the condition is also clear. 

In order to perform single electron pumping starting 
from a switch configuration opposite to the one discussed 
above, one has two similar conditions on Ti and 
Therefore, the necessary and sufficient condition for sin- 
gle electron pumping is eo < 2 minpTi, r 2 , T3, T4}. 

V. QUANTUM INTERFERENCE OF MFS 

The parity flipping in the NOT gate discussed above 
can be understood generally as a quantum interference 
of MFs. We consider a loop of 2N TSs with a vortex at 
the center of the common superconducting substrate. At 
the initial state, the odd/even constriction junctions are 
on/off, with no edge MF in the system. First, we turn 
off Ai, which produces two edge MFs 71 and 72 on TS(1) 
and TS(2) respectively. By turning on A 2 and A 2 at, and 
then turning off A3 and A2JV-I) and so on so forth, the 
MF 72 and 71 are transported to TS(21) and TS(21+1) 
respectively. Finally, we turn on X 2 i, which annihilates 
the two MFs. 



After the above sequence of switchings, the odd/even 
constriction junctions become off/on, opposite to the ini- 
tial configuration. Let us check the ground-state elec- 
tronic parity for the two configurations, keeping in mind 
that all the coupling constants Tj are positive when 
switched on as the topological property revealed ex- 
plicitly for the NOT gate with four TSs. The ground 
state \G) at the initial configuration is a fully occupied 
state by electrons defined by = (72.7-1 + «72j)/2 for 
1 < 3 < N. The electronic parity of the whole sys- 
tem P = (_i)£i=i.* 4* = n j= i,jvH72j72j-i) has the 
eigen value (— 1)^. 

For the final configuration, the ground state \G) is a 
fully occupied state by electrons defined alternatively by 
dj = (72j +«72j+i)/2 with 72W+1 = 71 ■ Its electronic 
parity can be calculated as 

P\G) = 1] (-^+i72i)|G) = (-1) N+1 \G), (5) 

3 = 1, N 

which is opposite to the one of the initial configuration. 
It is noticed that when the system of 2N TSs is iso- 
lated, the parity should be preserved, and thus the true 
ground state with reserved parity cannot be achieved. 
The attachment of a MF qubit or a QD assists the parity 
flipping as discussed above. 

It is easy to prove that without the central vortex, or 
generally with an even number of vortices, a system of 2N 
TSs takes the same ground-state parity after switching. 
Therefore, with an odd number of central vortices, the 
quantum interference of two MFs is constructive, which 
flips the ground-state parity, whereas with an even num- 
ber of central vortices, the quantum interference is de- 
structive leaving no parity change in the system. Quan- 
tum inference of MFs was discussed for a system with 
superconductor and two ferromagnets on the surface of 
topological insulator! 2 lj. The two phenomena share the 
same topological origin induced by the central vortex. 
However, in the present case, the parity is reversed upon 
the switching process, whereas in the case addressed be- 
fore an injected electron at one lead induces either emis- 
sion or absorbtion of an electron (strictly speaking, odd 
number of electrons), which preserves the parity. 

VI. DISCUSSIONS 

The mechanism of single electron pumping in the 
present proposal is different from previous ones in lit- 
erature [6|-|8(. Unlike the case where electrons are driven 
by biased potentials between QD and electrodes @, we 
switch the coupling between QD and the charge reser- 
voir, i.e. the NOT gate, instead of the bias. Since the 
present single electron pumping is independent of the de- 
tails of the switching process as far as it is slow enough, 
it is also different from the previous adiabatic pumping 
schemes where the details of pumping parameters dur- 
ing the operations have a significant effect on the direc- 
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tion and magnitude of the pumping current 0,9- In this 
sense, the present scheme can be called topological charge 
pumping. 

In our schemes, the manipulation of edge MFs is topo- 
logically protected by edge excitation gap which sets the 
limitation of the working temperature. For a sample 
size of 150xl50nm 2 , the edge excitation gap is estimated 
as ~0.01Ao with Ao the superconduction pair potential 
27]. For Ao=lmeV, the gap is ~ \QQniK, which is not 
hard to achieve experimentally in these days. It is wor- 
thy to note that there is a work by Akhmerov showing 
that even the occupation of higher in-gap excitations will 
not violate the topological protection [311 ] . 



duced by a vortex can be used to implement a quantum 
NOT gate for Majorana qubit. The working mechanism 
for the NOT gate can be understood as the Aharonov- 
Bohm interference of two Majorana fermions. Based on 
this phenomenon, we formulate a scheme for topological 
single electron pumping. Useful applications of these de- 
vices in quantum transport and quantum computation 
are expected. 



VII. SUMMARY 

In conclusion, we reveal that the mobility of Majorana 
fermion at edge of nano topological superconductor in- 
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